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In this article, for a radiative transport equation, we study inverse problems of 
determining a time independent scattering coefficient or a total attenuation by 
boundary data on the complementary sub-boundary after making one time input 
of a pair of a positive initial value and boundary data on a suitable sub-boundary. 
The main results are Lipschitz stability estimates. We can also prove the reverse 
inequality, which means that our estimates for the inverse problems are the best 
possible. The proof is based on a Carleman estimate. 



1 Radiative Transport Equation and Main Results 

We consider wave or particles propagating in a random medium. Let Q he a bounded 
domain of M'^, d > 2 with C^-boundary dQ. The scalar product in M'^ is denoted by 
a dot (•). Let V = V^,. = ( gf-j gf- ) • We let u{x,v,t) G M denote the angular 
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density at time t > and position x G M"^ with the velocity v E V, where V = 
{v eR'^;0<vo< \v\ < vi}. 

Let aa{x,v) and as{x,v) denote the absorption and scattering coefficients, respec- 
tively. Note that aa and as are positive measurable functions: 

: X 1/ ^ M, a^: ^xV ^R. (1.1) 

We introduce the total attenuation as at = aa + Cs- The following radiative transport 
equation, which is a linearized Boltzmann equation, governs u{x,v,t) for x E Q, v E V, 
<t <T, 

Pu := Pou{x,v,t) + at{x,v)u — as{x,v) / p{x,v,v')u{x,v\t)dv' = 0, (1.2) 

Jv 

where 

Pqu := dtu{x, v,t) + V ■ Vm(x, v, t). (1.3) 
The phase function p{x, v, v') satisfies 

p{x,v,v')dv' = 1 forall(x,f). (1.4) 

Equation (1.2) describes transport in a random medium such as light in biological tissue 
[H [2] , neutrons in a reactor [8], [TT] , and light in the interstellar medium [9] and atmo- 
spheres [29]. We let z/(x) be the outward normal unit vector to dfl at x G dfl. We define 
F-i- as 

T+ = {{x,v) EdnxV;u{x)-v>0}, = {{x,v) e dn x V; iy{x) ■ v < 0} . (1.5) 

We impose the following boundary conditions. 

u{x,v,0) = a{x,v), xeQ, V E V, (1.6) 
u{x,v,t) = g{x,v,t), 0<t<T, {x,v)eT^. (1.7) 

We consider inverse problems of determining at or ag by boundary data u{x,v,t), 
{x,v) G r+, < t < T after setting up an initial value (1.6) and boundary value (1.7) 
once. Our inverse problem is motivated by optical tomography, in which we recover 
at and as from boundary measurements (e.g., [HE])- An incident laser beam g{x,v,t) 
enters the sample on the boundary r_, and the outgoing light u{x, v, t) is measured on 
the boundary r+ x (0,r). 
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We refer to works concerning inverse problems on transport equation. ChouUi and 
Stefanov [10] proved the uniqueness of at and as for the scattering operator. We define 
the albedo operator A as 

A[g] = u{x, V, t), (x, v) eT+, <t <T, (1.8) 

by assuming that the initial value is zero. Stability in determining some coefficients 
among at, a^, p is proved by the angularly averaged albedo operator [1] and by the 
full albedo operator [5]. For the inverse problems in Refs. |1] and |5], the input-output 
operation can be limited to the boundary and the initial value can be zero, but one has 
to make infinitely many measurements. For the stationary transport equation, the non- 
uniqueness in the coefficient inverse problem with the albedo operator was characterized 
by gauge equivalent pairs in [31] , and the Lipschitz stability for gauge equivalent classes 
was proved for the time- independent radiative transport equation in [27]. See also review 
articles [31 ED] for coefficient inverse problems for the radiative transport equation. 

Klibanov and Pamyatnykh [23] proved the uniqueness of at by the boundary values 
of u. The formulation in [22] is different from [1], [S], [TU] and meaures a single output 
on r+ X (0, T) after choosing initial value and boundary data on r_ x (0, T). 

In this article, we adopt the same formulation as in [23] and we consider the inverse 
problems of detemining as or at by boundary value on r+ x (0,T) with a suitable 
single input of the initial value. Our main results are Lipschitz stability estimates in 
determining ag or at- To the best knowledge of the authors, there are no publications 
on the Lipschitz stability with a single measurement data. Bukhgeim and Klibanov [6] 
proposed a methodology for proving the uniqueness and the stability for inverse problems 
with a single measurement, on which [23] is based. Their method uses an L^-weighted 
estimate called a Carleman estimate for solutions to the differential equation under 
consideration. The Carleman estimate dates back to Carleman [7]. See Hormander [13] . 
Isakov [18], and Lavrent'ev, Romanov, and Shishat-skii [26]. As for inverse problems 
by Carleman estimate, we refer for example to Imanuvilov and Yamamoto [TS], [TB] , 
Isakov [T7], [19], Klibanov [20], [21], Klibanov and Timonov [21], and Yamamoto [32] . 
Moreover see Klibanov and Pamyatnykh [22] for the Carleman estimate for a transport 
equation and an application to the unique continuation, and Klibanov and Yamamoto 
|25j for the exact controllability for the transport equation. Prilepkov and Ivankov [2B] 
discusses an inverse problem of determining a t-function in the case where at depends 
on X, V, t. 
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Throughout this article, denotes usual Sobolev spaces. We set 

X = H\0, T- L'^in X V)) n H\0, T- L^{VL x V)). 

For arbitrarily fixed constant M > 0, we set 

U = {ueX- \\u\\x + \\Vu\\mio,T;L^{iixv)) < M}. (1.9) 

Now we are ready to state our main results. 
Theorem 1 (determination of at) 

Let u'^ — u{at){x,v,t), k — 1,2 be solutions to the transport equation: 

dtu{x, v,t) + V ■ Vu + cr^(x, v)u — as{x, v) / p{x, v, v')u{x, v', t)dv' = 0, 

Jv 

u{x, V, 0) = a{x, v), X & D,, V & V, k — 1,2, 
u^g onr_x(0,T). 

Let u'' e U and ||cr^||L°°(Qxi/)7 ||c'"s||L°°(fixy) < M. We assume that 

T>—sup\x-xo\ (1.10) 

with some Xq ^ fl and 



a{x,v)>0, {x,v)enxV. (1.11) 
Then there exists a constant C — C{M) > such that 
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/ {i^{x)-v)\dt{u'-u^){x,v,t)\''dSdvdty <\\al-aULHnxV) 



Theorem 2 (determination of ag) 

Let u'^ = u{ag){x,v,t), k = 1,2, be the solution to the transport equation: 

dtu{x,v,t) + V ■'S7u + at{x,v)u — ag{x,v) / p{x,v,v')u{x,v',t)dv' — 0, 

Jv 

u{x, V, 0) = a{x, v), X & fl, V & V, 
u = g on r_x(0,r), A; = 1,2. 

Let & U and \\(Jt\\L°°{Q.y.v), ||L°°(nxy); k = 1,2. We assume (1.10) and (1.11). 
Then there exists a constant C — C{M) > such that 

/ Hx)-v)\dt{u^-u^){x,v,t)\''dSdvdty\\\al-al\\L^n^v) 



In (1.12) and (1-13), the second inequalities show the Lipschitz stabihty for the inverse 
problems, while the first inequalities are related to the initial/boundary value problems 
in whcih we are required to find dtu on r+ x (0,T) for given a and {a^,as), {at,ag), 
k = 1,2. We obtain both-sided estimates and so the estimates for the inverse problems 
are the best possible. 

For the Lipschitz stability for the inverse problems, we need the positivity (1.11) up 
to the boundary d{Q x V) of the initial value. Measurements must be set up so that this 
positivity is guaranteed. The posivitiy condition is restricting but can be achieved in 
practice for example as follows. Let us consider optical tomography of the human brain 
(cf. [121 E]). We use a continuous- wave near-infrared laser beam and modulate the light 
by using an optical device. Before being temporally varied, the time-independent light 
is applied to the head. The light is then scattered in different directions in the brain, 
and comes out. Thus, in this setup, we can consider that the initial angular density 
a{x, v) in the head is positive in Q x V. 

Moreover we have to assume (1.10), that is, the observation time T should be large 
compared to the size of the domain Q. This is a natural condition because the transport 
equation has a finite propagation speed, which can be seen by (1.3). 

In order to prove Theorems 1 and 2, it is sufficient to prove the linearized inverse 
problem. More precisely. 
Theorem 3 

We consider 

dtU+v-Vu+atU—Gs / p{x,v,v')u{x,v' ,t)dv' = f{x,v)R{x,v,t), x & il, v G V, < t < T, 
Jv 

u{x, v,0) = 0, X G Q, V & V, 

u = on r_x(o,r). 

We assume 

R, dtR G L\0, T; L°°{n x V)), at, a, G L°°(^] x 1/), 

and 

u G //^(O, T; L^iVl X V)) n i/^(0, T; L^Q x V)), Vu G //^(O, T; L^iVt x V)). 
We further assume (1.10) and (1.11). Moreover let 

R{x,v,Q)>{), {x,v)eTTxV. 
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Then there exists a constant C > such that 

(1.14) 

for any f G x V). 

In fact, for the proof of Theorem 1, setting u = — v?, f = a] — and R = —u^, 
we have the above hnearized inverse problem. By the regularity assumption of u^,u'^, we 
can apply Theorem 3 to obtain the conclusion (1.12). We can similarly derive Theorem 
2 from Theorem 3. 

The article is composed of 4 sections. In section 2, we prove the first inequality of 
(1.14). In section 3, we prove a key Carleman estimate and in section 4, we complete 
the proof of Theorem 3. 



2 Proof of the first inequality in Theorem 3 

We set 

Q^nxV. 

In this section, C > denotes generic constants which are independent of /. 
We prove the first inequality by energy estimation. We set 

Ui — dtu. 

Then 

Pui{x, V, = f{x, v)dtR{x, V, C), {x,v) e Q, < i < t. 

Multiplying Pui{x, v, ^) — fdfR by 2ui{x, v, ^) and integrating over (x, v,^) & Qx (0, t), 
we have 

f dt( I \dtu{x,v,i)\'^dxdv\di+ I [ v-V{\dtu{x,v,^)\'^)dxdvd^+ [ [ 2at\dtu\'^dxdvd^ 
Jo \Jq J Jo Jq Jo Jq 

(2.1) 



\JQ 



2asdtu{x,v,^) p{x,v,v')dtu{x,v' ,^)dv'^ dxdv^ d^ 

+ / / 2{dtu){x,v,^)f{x,v){dtR){x,v,^)dxdvd^. 
Jo Jq 
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Integrating by parts and noting that at G L°°{Q), for < t < T, we have 

vV{\dtu\^)dxdvd^^ III {vv)\dtu\^dSdvdi 
I Jo Jv J on 

{vv)\dtu\^dSdvdi> [ [ {vv)\dtu{x,v,i)\^dSdvdi 

Jo Jv^ 



\^r+ JT-/ Jo Jr^ 

and 



2\at\\dtu\^dxdvdi<C I I \dtu{x,v,C)?dxdvd^. 

Jo JQ 

Since u{x, v, ^) = dtu{x, v, r])dr) by u{x, v, 0) — 0, we have 



I [the first term on right-hand side of (2.1) 
< C 







\dtu{x,v,^)\ i^j \dtu{x^v\^)\dv'^ dxdv^ d^ 



dtu{x, v', r])dr] 



dv' \ dxdv ] d^ 



= c f ( [ \dMx,v,o\( [ 

Jo \JQ \JV 

< ^ ( / / \dtu{x,v' ,r])\drjdv' ] dxdvd^ 

Jo Jq \Jv Jo / 

< c[ [([ [\\dtu{x,v,0\^ + \dtu{x,v',rj)\^)dridv']dxdvd^ 
Jo Jq \Jv Jo J 

[ f \dtu{x,v,i)\^dxdvdi, 0<t<T. 

Jo Jq 



iQ 

< C 



'0 JQ 

At the second last inequahty, we used 



\dtu{x,v,^)\\dtu{x,v',Ti)\ < \dtu{x,v,i)\^ +\dtu{x,v',Ti)\'^ 



Therefore, since 



' dt i \dtu{x,v,^)\'^dxdv \ d^ = / \dtu{x,v,t)fdxdv— / \f{x,v)R{x,v,Q)\'^dxdv 
\Jq / Jq Jq 

by u{x, V, 0) = 0, applying the Cauchy-Schwarz inequality to the second term on the 
right-hand side of (2.1), we obtain 

/ \dtu{x,v,t)\'^dxdv (2.2) 
Jq 



< 



\dtu{x,v,t)\'^dxdv + / / {u ■ v)\dtu{x,v,^)\'^dSdvd^ 
Jo Jr+ 



< c[ [ \dtu{x,v,0\^dxdvd^ + C [ I \f\^\dtR\^dxdvd^ 

Jo JQ JO JQ 

+ \ f{x,v)R{x,v, 0)1'^ dxdv, 0<t<T. 

JQ 

By i?(-,-,0) e L'^iQ) and dtR G L'^{0,T; L°°{Q)), applying the Gronwall inequality to 
(2.2), we obtain 

/ \dtu{x,v,t)\dxdv<C\\f\\l,^Q^, 0<t<T. (2.3) 
JQ 

Applying this to the right-hand side of (2.2), we obtain 

i^ix) ■ v)\dtu{x,v,0\'dSdvd^ < 0<t<T. 



10 Jr^ 

Thus the proof of the first inequality is completed. 

3 Key Carleman estimate 

We set y — y{x, t), u — u{x, v, t) and 

Ly{x,t)^dty + v-Vy, x e t > 0, (3.1) 

where v & V is arbitrarily fixed. We define 

(p{x,t) ^\x-xo\'^ - (3.2) 

where < (3 < Vq. Then we can prove the key Carleman estimate. 
Lemma 3.1 

(i) There exist constants sq > and C > such that 

s\y{x,t)\''e^'^dxdt 



T Jn 

< c[ [ \Ly\'^e'^"^^'''^Uxdt + Cs [ [ {vv)y'^e^"^^''^^USdt 

J-T JQ. J-T Jvu(x)>0 

for all s > So and y e L'^{-T,T; H\n)) n H^{-T,T; L'^{n)) such that y{x,t) = 0, 

0<t<TifxEdQ satisfies ^{x) ■ v <0. 

(a) There exist constants Sq > and C > such that 

[ s\u{x,v,t)\^e'^"^dxdvdt 

-T JQ 



< C [ [ \Lu\^e^''^^''''^dxdvdt + Cs [ [ {u-v)u^e^''^^''''^dSdvdt 
J-T Jq J-t Jr+ 

for alls > So andy e L^{-T,T; L\Q))r\H^{-T,T; L\Q)) such thatVy G L\-T,T; L^{Q)) 
and y = on r_ X (— T, T) . 

Proof The part (ii) is proved by integrating the conclusion of (i) and so it suffices to 
prove (i). We set z{x,t) = e"^(^'*)y(x, t) and (L,z)(x,t) = e"^(^'*)L(e-*'^2). Then 

LsZ = dtz + V ■ Vz — s{{dt(p) + {v ■ V(p))z := dtz + v ■ Vz — sA{x, t)z. 

Here we set A = dtf + [v ■ Vf) = — 2/3t + 2v ■ {x — xo). Hence by y\r_ = 0, we have 

[ [ \Ly\^e^'''^''''Uxdt= [ [ \Lsz\^dxdt 
J-T Jn J-T Jn 

> -2s [ [ A{dtz + V ■ Vz)zdxdt = s [ [ {Adtiz"^) + Av ■ V{z^))dxdt 

rT r rT 



: [ [ {dtA + VA - v)z'^dxdt - s [ [ A{p-v)z^dSdt 
J-T Jn J-T Jan 



i-T f pT 
1 2 n\ I / ^2 



2s{\v\'-l3) / / z^dxdt-s / / A{u-v)z^dSdt 
J-T Jn J-T Jan 

> 2s{v^-/3) [ [ z^dxdt-Cs [ [ {yv)z^dSdt 

J-TJn J -T J {u{x)-v)>0 

Substituting z = e^'^y, we complete the proof. 



4 Proof of Theorem 3 

The proof is similar to Imanuvilov and Yamamoto [T3], [TB]. We set 

Ml = dfU. 

Then we have 

Pui = f{x,v)dtR{x,v,t), {x,v) e Q, < t < T, 
ui{x,v,0) = f{x,v)R{x,v,0), {x,v) G Q. 



We extend y,as,crt,p as follows. 

y{x,v,t) ■ 



dtu{x, v,t), t > 0, 
dtu{x, —V, —t), t < 0, 
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, , , (Js{x,v), t>0, \ at{x,v), t>0, 

as{x,v,t) = { at{x,v,t) = < 

-as{x,-v), t<0, y -at{x,-v), t<0, 

and 

p{x, V, v'), t > 0, 



p{x,v,v',t) 



—p{x, —V, v'), t <0. 



We set R = dfR, and we extend R by R{x,v,t) = —R{x, —v, —t) for t < 0. Therefore, 
for t < 0, we have dty{x,v,t) = —{d^u){x,—v,—t) and Vy{x,v,t) = V dtu{x , —v , —t) , 
and so 

y e H\-T,T;L'{Q))nL'{-T,T;L'{Q)), Vy G L\-T,T; L'{Q)). 
Moreover 

dty + V -Vy + aty - as pydv = f{x,v)R{x,v,t), {x,v) e Q, -T < t < T. (4.1) 

Jv 

By (1.10) we can choose ^ > such that ;^ < 1 and 

sup^gn \x-xo\ 

^ ■ (4.2) 

Later, we will bring 4- to 1 as close as possible. 
Therefore we have 

ip{x,±T)^\x-xo\'^ - ^T^ <0, xen 

and 

ip{x,o) > 0, X eil 

by Xo ^ Q. Hence we can choose S > sufficiently small such that 
ip{x, t) <-d, -T < t < -T + 25, T - 25 < t < T, X 

and 

cp{x, t) > 5, -5 < t < 5, X e ^. 
We fix X e C^(M) such that < x < 1 and 



x{t) 



1, -T + 25<t<T-25, 

0, -T < t < -T + 5, T - 5 < t < T. 
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We set 

Then we have 

dtz + V ■ V z + atz — as / p{x,v,v' ,t)z{x,v' ,t)dv' (4.3) 

Jv 

= X{t)f{x, v)R{x, V, t) - {dtx)y, {x, v)eQ,-T<t<T, 

Z{-,;±T)^0 

and 

dfW + V ■ Vw + atw — Gs / p{x,v^v' ,t)w{x,v' ^t)dv' (4.4) 

Jv 

= Xe'^R - {dtx)ye''' - s{dtip + {v ■ V^))w, {x, v)eQ,-T<t<T, 

W{;;±T)=0. 

Henceforth C > denotes generic constants which are independent of s > and /. 

By at e L°°{Q X (-T, T)) and p e L'^ifl x V x V x (-T, T)), we see that 




TJQ 



as{x,v,t) / p{x,v,v' ,t)z{x,v' ,t)dv' 



V 



e'^"^^'''^Uxdvdt (4.5) 



< I Ws{x,v,t)\^\p{x,v,v',t)\^\z{x,v',t)\^dv'^ e^''^^'-^''^ 

< ^ J J (^J (^J Hx,v',t)\^dv'^d^e^"^^''^^^dxdt<C j J \z{x,v,t)\^e^''^^'''^'^dxdvdt 



and 



\at{x,v,t)z{x,v, t)|V^'^(^'*)(ixM< C / / \z{x,v,t)\^e'^"^^''''^dxdvdt. 

J-T Jo 



-TJQ J-TJQ 

Therefore, applying Lemma 3.1 (ii), we have 



S 



+ c ' 



j j \z\''e''^^''''^dxdvdt<C j \f{x,v)\^ (^j e-^'^''d?je''''^^=''''Uxdv 

\dtx\Ye^'^^'''^^dxdvdt + Ce^' j j {v ■ v)y^dSdvdt 

J — T J r_i_ 



-TJQ 
"T 

+ C I I \z\^e^''^^''^^^dxdvdt 




TJQ 
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for all large s > 0. By absorbing the last term on the right-hand side into the left-hand 
side by choosing sufficiently large s > 0, and we obtain 



s J \z\''e^"^^'''''^dxdvdt <C J \f{x,v)\^ (^j' e-^'^''d?j e'^'^^'^'^'Uxdv (4.6) 

+C [ [ \dtxWe^"^^''''^dxdvdt+Ce^' [ I {u ■ v)y^dSdvdt. 
J-T Jq J-t Jr+ 

Multiplying (4.4) by —2w{x,v,t) and integrating over Q x (0,T), we obtain 



T /■ rT 

2 



dt{w )dxdvdt — / / v-\/{w )dxdvdt 

JQ Jo 



— / / 2atW^dxdvdt — 2as / / I / pw{x,v' ,t)'w{x,v,t)dv' \ dxdvdt 
Jo Jq Jo Jq \Jv J 

If 2{xe"^fRw - {dtx)ye"^w - s{dt(p + v ■ Wip)w'^)dxdvdt. 
Jo Jq 

w — x{dtu)e^'^ for f > 0, we have ^^; = on r_ x (— T", T). Therefore we have 



[the left-hand side] 



/ w{x,v,QiYdxdv — III w'^{v ■ v)dSdvdt 
Jq Jo JdQ Jv 

J w{x,v,Oydxdv - j (^j +j ^w^{v ■v)dSdvdt 



Qfdxdv- I I w^{iy ■v)dSdvdt 
Jq Jo Jr+ 

{x,v,0)\^e^"^^'''°Uxdv- [ [ x'^\dtuWiy-v)e'^''^^'''^^dSdvdt. 

Jo Jr+ 



\y[ 

By the Cauchy-Schwarz inequality and an argument similar to (4.5), we have 
I [the right-hand side] | 

cT r rT 



< Cs [ [ z^e'''^^''''^dxdvdt+C f f z''e'''^^'=''^dxdvdt 
J-T Jq J-t Jq 

+ C [ [ {dtx)V(^^"^^''''^dxdvdt + C [ [ fe^"^^''''Uxdvdt. 
J-T Jq J-t Jq 



-TJQ J-TJQ 

Therefore 

rT r rT 



\dtu{x,v,0)\^e^"^^''''^^dxdv <C [ [ \f\^e^"^^''''^dxdvdt+C [ [ {dtx)Ye^"^^''''^dxdvdt 

J-T Jq J-t Jq 
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+Cs [ [ z^e^''''^^''Uxdvdt+ [ [ \dtuWu ■ v)e^''''-''''USdvdt. (4.7) 

J-T Jq J-T Jt+ 

Substituting (4.6) into the third term on the right-hand side of (4.7), we obtain 
Jq 

< cj \fix,v)f(^j e-^"^''d?je^"^^'''''Uxdv + C J' J \dtxWe^'^^''''^dxdvdt 
+ Ce^' I I {v v)y^dSdvdt. 

Here by the definition of X) we see that dtX 7^ only if —T + 8 < t < —T + 25 or 
T — 25<t<T — 5, which imphes (p[x, t) < —S. Hence 



\dtx\%^e^''^^''''Uxdvdt < Ce-^'^ I [ y^dxdvdt. 

J-T Jq 



-TJQ J-TJQ 

By the Lebesgue theorem, we have 



j e-^'^^'dt = o(l) as s ^ oo. 



By ■u(-,-,0) = 0, we obtain dtu{x,v,0) — f{x,v)R{x,v,0). Consequently, R{x,v,0) ^ 
for {x, v) e Q, we obtain 

/ \f{x,v)\^e^"^^'''^^dxdv (4.8) 
Jq 



< o(l) / \f{x,v)\^e^'^^^'''^dxdv + Ce-'''^ [ I y^dxdvdt 
Jq J-t Jq 



'Q 

+ Ce^' I I {v v)y^dSdvdt 
J — T J r_i_ 



for all large s > 0. By (2.3), noting that y is the extension of dtu to t < 0, we have 

/ \y{x,v,t)\^dxdv < C||/|U2(Q), -T<t<T. 
Jq 

Note that \x — xo\ > s > 0, x e Q with some constant £ > because xq O,. By 
(4.8) we obtain 

{l-o{l))e^'' [ \f{x,v)\^dxdv <{l-o{l)) [ |/(x,v)|V*'^(^'°)cixciw 
Jq Jq 

< Ce-^'^T [ \ f{x,v)\^dxdv + Ce^' [ [ {u- v)y^dSdvdt 
Jq J-t Jr+ 
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for all large s > 0. Since y is the extension of dtu to t < 0, we have 

{(1 - o{l))e^'' - Ce-'^'^T} I \ f{x,v)\^dxdv < Ce^' I [ {u ■ v)\dtu\^dSdvdt 

JQ Jt+ Jo 

for all large s > 0. Choosing sufficiently large s > 0, we complete the proof. 
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